NASA TN D-7538 



NASA TECHNICAL NOTE 


CO 

un 

i 


■C 



.'fNASA-TN-D-7538) IBIAMGULATIGM OF 

moliistatiom cabeba data to locate a 
COBVED LIME IH SPACE {MAS A) 42 p BC 
;$3.25 CSCL 


14E 


Hi/30 


y 74-22478 


Onclas 

38084 



TRIANGULATION OF 
MULTISTATION CAMERA DATA 
TO LOCATE A CURVED LINE IN SPACE 


by Cliffofd L. Fricke 

Langley Research Center 
Hampton , Va. 23665 



NATIONAL AERONAUTICS AND SPACE ADMINISTRATION • WASHINGTON, D. C. ♦ MAY 1974 



1. Report No, 2. Government Accession No. 

NASA TN D-7538 

4 . Title and Subtitle 

TRIANGULATION OF MULTISTATION CAMERA DATA TO 
LOCATE A CURVED LINE IN SPACE 

7. Author(s) 

Clifford L. Fricke 


9. Performing Organization Name and Address 

NASA Langley Research Center 
Hampton, Va, 23665 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546 

15. Supplementary Notes 


3. Recipient's Catalog No. 


5. Report Date 


6. Performing Organization Code 

May 1974 


8. Performing Organization Report No. 

L-9268 


10. Work Unit No. 

879 - 11 - 36-01 


11. Contract or Grant No. 


13. Type of Report and Period Covered 

Technical Note 

14. Sponsoring Agency Code 


16. Abstract 


A method is described for finding the location of a curved line in space from local 
azimuth as a function of elevation data obtained at several observation sites. A least- 
squares criterion is used to insure the best fit to the data. The method is applicable to the 
triangulation of an object having no identifiable structural features, provided its width is 
very small compared with its length so as to approximate a line in space. The method was 
implemented with a digital computer program and was successfully applied to data obtained 
from photographs of a barium ion cloud which traced out the Earth's magnetic field line at 
very high altitudes. 


17. Key Words (Suggested by Author (s)) 

Triangulation, multistation 
Ion cloud 

Stereotriangulation 
Photogrammetry 
Ballistic camera 
Multistation triangulation 


19. Security Cl a ssif. (of this report) 
Unclassified 


18. Distribution Statement 

Unclassified - Unlimited 


20. Security Classify (of this page} 

21. No. of Pages 

Unclassified 

*3 | 


STAR Category 30 


22. Price* 

$3.25 


For sale by the National Technical Information Service, Springfield, Virginia 22151 

t 

! 
















TRIANGULATION OF MULTISTATION CAMERA DATA TO LOCATE 
A CURVED LINE IN SPACE 

By Clifford L. Fricke 
Langley Research Center 

SUMMARY 

A method is described for finding the location of a curved line in space from local 
azimuth as a function of elevation data obtained at several observation sites. A least- 
squares criterion is used to insure the best fit to the data. The method is applicable to 
the triangulation of an object having no identifiable structural features, provided its 
width is very small compared with its length so as to approximate a line in space. The 
method was implemented with a digital computer program and was successfully applied 
to data obtained from photographs of a barium ion cloud which traced out the Earth's 
magnetic field line at very high altitudes. 

INTRODUCTION 

A unique and powerful tool in magnetospheric studies involves the deposition of 
barium vapor at a point in the magnetosphere through the use of chemicals or explosives 
carried aloft by rockets. The barium atoms are rapidly ionized by sunlight and thus form 
a barium ion cloud which extends along the magnetic field line and becomes "frozen" to it. 
The barium ion cloud, on account of its resonant scattering of sunlight, is visible to ground 
sites when viewed against the night sky and hence serves to delineate the magnetic field 
line over a considerable arc length. This condition permits a determination of magnetic 
field line orientation and shape. From the motion of the cloud one may obtain the convec- 
tive motion of magnetospheric plasma and hence the electric fields which drive such 
motions. 

The accuracy of locating an object in the distant magnetosphere by triangulation 
methods is extremely limited by the relatively short baselines available. In the case of 
the barium ion cloud (ref. I), even though the observation sites were widely dispersed on 
the Earth, the cloud altitude was about five times the baseline distance. Hence, it was 
necessary to use every means possible to improve the accuracy such as calibration for 
distortion in the cameras, orientation of cameras using stars in the photographs, and the 
use of more than two observation sites. 
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Existing triangulation methods were limited in one respect or another: methods 
using several observation sites triangulated only on points (refs. 2 and 3) and methods 
which triangulated on lines could use only two observation sites (refs. 4 to 6). Table I 
shows a comparison of the various methods. 

In view of these limitations, it was desired to develop a triangulation method which 
allows for a cloud that is curved in space, extends over a large arc, and is to be photo- 
graphed at several observation sites. The difficulty in obtaining a solution by triangu- 
lating on an extended object from several sites lies in the fact that the intersection of the 
several surfaces defined by the necessarily inaccurate data from several sites does not 
determine a unique line in space. Thus, some criterion is needed to determine the most 
probable solution when all the data are taken into account. 

Methods already exist (refs. 7 to 9) for precise conversion of images on photographs 
to pointing directions from each site, such as azimuth as a function of elevation. This 
note describes a method and a computer program that find the most probable line solution 
from input data on azimuth as a function of elevation from several observation sites. 

FORTRAN VARIABLES AND SYMBOLS 

FORTRAN variables are the same as algebraic symbols, except that they are under- 
lined when used in equations in the text. 

A equatorial radius, 6378.166 km; in subroutine SUMRES, angle in radians used 

to weight residuals 

AZ(L,N) azimuth angle from site L, data point N 

a,b,c coefficients in quadratic equation 

B polar radius, 6356.784 km 

BC(I,L,N) Ith coefficient in least-squares fit for site L, data point N 
CLAT geocentric latitude 

CR geocentric radius 

C1,C2,C3 coefficients in least-squares fit to data for azimuth as a function of elevation 
(az-el) 
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D 

D(L) 

DA 

DC 

DDA 

DEG 

DIST(L) 

DNA 

DR 

DX,DY,DZ 

E 

E(N) 

EL(L,N) 

EM 

E1,E2,E3 

F 

GLAT 


residual (minimum angle between point P and az-el curve) 
residual from site L 
increment given to PLON 
discriminant of cubic equation 

interpolation correction given to PLON to give minimum point 
degrees per radian 

angle between point P and end data point on az-el curve, site L 
altitude increment 
increment given to PLAT 

axes in topocentric coordinate system; on horizontal plane in easterly direc- 
tion, on horizontal plane in northerly direction, and in vertical direction 
perpendicular to horizontal plane, respectively 

root- mean- square residual; in subroutine GGRGCN, ratio A 2 /B 2 

minimum root-mean-square residual with respect to both latitude and longi- 
tude; thus, the solution at point N is defined 

elevation angle from site L, data point N 

minimum root-mean-square residual (with variable longitude) for a given 
latitude 


in subroutine LONMIN, root-mean-square residual for three consecutive 
longitudes 

flattening factor, 1/298.3; in subroutine SUMRES, weighting factor depending 
on DIST(L) 

geographic latitude 
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GLON 


H 

HI 

IN 

IS 

L 

LP 

N 

NBL 

NB(L) 

NC(L,N) 

ND 

NO 

NP(L) 

NS 

NV 

N1 

N2 

P,Q,R 


geographic or geocentric longitude 

geographic altitude, km 

altitude of first estimated trial point 

number of iterations in subroutine MINISOL (usually 3) 

number of sites using data points on the end of az-el curve 

observation site number 

index to designate line (LP=1) or point (LP=2) solution 

in program LARC, data point number, or solution point number; in subroutine 
RESDUE, data point nearest the trial point (PAZ,PEL) 

number of az-el data points (used for NB(L)) 

number of az-el data points from site L 

az-el data point number of site L used to obtain solution point N 

index used to increment solution point number N by ±1 

number of az-el data points used in least-squares fit 

data point number nearest solution point from site L 

number of observation sites 

particular data point within the set NO 

smallest solution point number 

largest solution point number 

coefficients in cubic equation 
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PAZ 

PEL 

PI 

PLAT 

PLATG 

PL ATI 

PLON 

PLON1 

PR 

R 

RA 

RAD 

RP 

SH(L) 

SLAT(L) 

SLON(L) 

SR 

TALT(N) 

TLAT(N) 

TLON(L) 


azimuth of point P 
elevation of point P 

TT 

geocentric latitude of point P 
geographic latitude of point P 

geographic latitude of first estimated trial point P 
longitude of point P 

longitude of first estimated trial point P 
geocentric radius of point P 
radius of revolution of spheroid 
range from observation site S to point P 
radians per degree 

projection of geocentric radius on equatorial plane 
altitude of observation site L 
geographic latitude of observation site L 
longitude of observation site L 
radius from Z-axis of observation site 
altitude of solution point N ' 

geographic latitude of solution point N 
longitude of solution point N 



WT(L) weighting factor for observation site L 

X, Y,Z geocentric coordinate axes; in equatorial plane in direction of Greenwich, in 

equatorial plane directed to make a right-handed system, and in direction 
of north pole perpendicular to equatorial plane, respectively 

XA,YA,ZA geocentric components of line-of-sight vector from observation site S to 
trial point P 

XD angle between point P and data point on az-el curve 

XP,YP,ZP geocentric components of trial point P 
XS,YS,ZS geocentric components of observation site S 

XI, Yl independent and dependent variables, respectively, in least- squares curve 

fit to az-el data 

0,0, 0 Eulerian rotation angles 

GENERAL DISCUSSION OF MULTISTATION TRIANGULATION 
Case of a Point Object 

Consider first the simple case of a point object in space observed from several 
sites. (Fig. 1 illustrates the case of three sites.) In general, the measured lines of 
sight will not intersect because of errors in measuring the pointing direction from each 
site. One would expect the probable error in pointing to be the same at each site, where 
the error in pointing is defined as the angle between the measured direction and the actual 
direction of the object. In order to find the most probable solution, consider a trial solu- 
tion point P in space and define the residual from a site as the angle between the mea- 
sured direction from that site and the direction of the trial solution (D(l), D(2), and D(3) 
in fig. 1). The most probable solution, if it is assumed that pointing errors at each site 
are random, then, would be the one which minimizes the sum of the squares of the resid- 
uals D(l) 2 + D(2) 2 + D(3) 2 . (See ref. 10, pp. 107-109.) 

Case of a Curved Object 

Next, consider the problem of locating a curved line object in space which is illus- 
trated in figure 2 for three observation sites. The curve for the azimuth as a function of 
elevation from a single site defines a conical surface in space. With data from only two 
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sites, intersection of the surfaces is unique, but with three or more sites, the intersec- 
tions are no longer unique. 

Generally, there is no distinguishable feature on any part of the cloud. The end 
points cannot be determined since they fade out gradually, and hence the location of the 
end depends on the exposure of the photograph, and brightness of the sky background. 

The center of the cloud length may sometimes be brighter, but it cannot be located accu- 
rately. Hence, if the curvature of the azimuth as a function of elevation is small, only 
errors perpendicular to the curve are important, since any error along the curve merely 
slides the curve on itself. It is reasonable to expect, then, that the probable error per- 
pendicular to the surface defined by the az-el curve is the same for each site. Thus the 
residual of a trial point will be defined as the angle between the line of sight of that point 
and its projection on the conical surface defined by the curve for the variation of azimuth 
with elevation from the specified site. This angle, of course, represents the minimum 
angle between the line of sight to the trial point and any line of sight on the conical sur- 
face (D(l) in fig. 2). If the trial point is to lie on the most probable solution curve, then 
the sum of the squares of the residuals from all sites must be minimized. (This pro- 
cedure is equivalent to minimization of the root-mean-square (rms) value of the 
residuals.) 

It has been assumed throughout that the object is a line in space. If the object has 
lateral dimensions, then it must possess a center line which is identifiable as such from 
every observation site for the method to be applicable. 

Minimization of Root-Mean- Square Residuals 

Even though the desired solution is a curved line in space, the method described in 
this paper successively solves for specific points on that line. For convenience, the inde- 
pendent variable is altitude (which will be converted to geocentric radius) so that each 
point will be at a selected altitude; if the cloud had extended beyond the magnetic equator, 
then the independent variable should be latitude in order to be single-valued. 

Thus, consider a trial solution point P which will remain at a fixed geocentric radius 
(which is essentially equivalent to a fixed altitude); the problem is to find the latitude and 
longitude which minimize the root- mean-square residuals. First, assume a trial latitude 
A and vary the longitude along A until a minimum root- mean- square residual EM^ is 
found, as indicated in figure 3. Next, increment the latitude to some new value B, and 
again vary longitude until a minimum as indicated by EMg is obtained. Thus EM, the 
minimum root- mean- square residual with variable longitude can be obtained as a function 
of latitude, and one can readily find the latitude that minimizes EM; this latitude and its 
corresponding minimizing longitude is the point at the selected altitude that lies on the 
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most probable solution curve. A description of the computer programs which use this 
method to find the most probable solution curve is given in the following sections. 


COMPUTER PROGRAMS 

The main program, called LARC (listing and flow chart are given in appendix A), 
in addition to reading cards, processing data, and printing results, has the main function 
of iterating the altitude at desired preselected intervals. Program LARC calls subroutine 
MINISOL which first calls subroutine LONMIN in order to minimize the root- mean-square 
residuals with respect to longitude, and then iterates the latitude in order to find the min- 
imizing latitude. 

The subroutines will be discussed first and are listed in appendix B. FORTRAN 
variables in the computer programs will have the same name in the text, but when used 
in algebraic equations, they will be underlined. 

Appendix C describes the geographic, geocentric, and topocentric coordinate sys- 
tems. The conversion from geographic to geocentric (subroutine GGRGCN), as well as 
the formulas for the inverse conversion (subroutine GCNGGR), are also given in 
appendix C. 


Conversion From Geocentric to Topocentric Coordinates 

Subroutine PAZEL converts a trial solution at a point P (given in geocentric coordi- 
nates, PLAT, PLON, and PR) to local azimuth and elevation from an observation site S 
with geographic coordinates SLAT, SLON, and SH. Figure 4 applies. 

The first part of the program converts the observation site to geocentric coordinates 
XS, YS, and ZS by using the same method as in subroutine GGRGCN. 

The geocentric coordinates of the trial point are converted to Cartesian coordinates 
XP, YP, and ZP. The geocentric coordinates of the line-of-sight vector from observation 
site S to trial point P are 


XA 


XP - XS 

YA 

= 

YP - YS 

ZA 


ZP - ZS 

J 


. 


The geocentric axes X, Y, and Z can be brought into coincidence with axes DX, DY, 
and DZ by rotation through the Eulerian angles 
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0 = SLON + 90° 


9 = 90° - SLAT 
0 = 0 

Thus, from use of formulas derived in reference 11, 


DX 


-sin SLON cos SLON 0 


XA 

DY 

- 

-sin SLAT cos SLON -sin SLAT sin SLON cos SLAT 


YA 

DZ 


cos SLAT cos SLON cos SLAT sin SLON sin SLAT 


ZA 

- 


- 


- 


Then the azimuth angle is 


DX 

AZ = tan' 1 -=- 
DY 


and the elevation angle is 


EL = tan -1 


DZ 


DX + DY 



Sorting of az-el Data 

Since certain subroutines require the az-el data to be ordered from one end of the 
curve to the other, it is necessary to insure that they are. Subroutine SORT(L) sorts the 
az-el data from each site L, for convenience, in such a way that the first point from each 
site corresponds to the high altitude end of the cloud; this is accomplished by setting the 
variable 


XS= -(-1) L 

This particular equation, of course, was made to hold for a particular orientation of the 
cloud and a particular set of observation sites. 
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Since elevation was the independent variable and single -valued, the ordering was 
done in terms of that variable only. 

Computation of Second-Order Least- Squares Fit to 
Successive az-el Curve Segments 

The accurate computation of the residual is an important part of the triangulation 
method. The original az-el data, in a striving for accuracy, usually consists of many 
arbitrarily located closely spaced points along the curve, so that a least-squares curve 
fit can be used to reduce the random errors involved in measuring the cloud center line. 

A second-order fit is sufficient for defining a segment of the curve since only a short 
interval is needed in the vicinity of the trial point P. The number of points used in the 
curve fit is NO and depends on the number and quality of data points. 

Subroutine BCOEF(L,N) computes the three coefficients BC(3,L,N) for a second- 
order least-squares curve fit from a given observation site numbered L, with the data 
point N as an origin. The NO data points used in the curve fit are centered about the data 
point N (except near the ends of the curve) and hence NO is an odd number. It is to be 
noted that the coefficients are calculated for every data point, of which there are NB(L). 

An exact coordinate conversion to the point N as origin would involve the Eulerian 
angles AZ and EL, but since the angular deviations from this origin will always be small, 
the two orthogonal angular components are 

XI = AEL 


and 


Y1 = A AZ cos ( EL ) 

as illustrated in figure 5. The independent variable in the quadratic formula for the par- 
ticular data point NV within the set NO is 

XI = EL(L,NV) - EL(L,N) 


The dependent variable is 


Y1 = (AZ(L,NV) - AZ(L,N))cos 


EL(L,N) + EL(L,N) 
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Computation of Residuals 


Subroutine RESDUE(PAZ,PEL,L,D,XD) finds the residual D of a trial solution point 
when the azimuth PAZ and the elevation PEL of the line of sight from site L to the trial 
point P are given. 

The first part of the subroutine finds the data point N which is nearest the trial 
point PAZ, PEL, starting with the point from the previous calculation which is stored in 
NP(L) as a first try. 

The next part of the subroutine computes the residual. Point P has the coordinates 
(with data point N as the origin) 


XI = PEL - EL(L,N) 


PEL + EL(L,N) 

Yl = ( PAZ - AZ(L,N))cos = — — 

2 

The second-order curve fit to the data points is 
Y = Cl + C2 X + C3 X 2 


where 


Cl = BC(1,L,N), . . . 

The distance between point P and the curve is given by 
D 2 = (XI - X) 2 + (Yl - Y) 2 

To find the minimum distance, differentiate and set the result equal to zero: 


dY 

(XI - X) + (Yl - Y) — = 0 
— “ “ dX 


Substitution for Y and dY/dX from the preceding equations gives the standard form 
for a cubic 


11 



where 


X 3 + PX 2 +QX + R = 0 


1.5C2 
P = — — 
C3 


1 + C2 2 + 2C3(C1 - Yl) 

2C3 2 

C2(C1 - Yl) - XI 
R = ~ - ~ — 

2C3 2 

The solution of this cubic equation is standard, but the type of solution depends on 
the value of the discriminant DC (as shown in the listing of subroutine RESDUE). If DC 
is greater than zero, there is one real root which is computed. If DC is less than zero, 
there are three real roots. The residual is thus 

D = ^(Y - Yl) 2 + (X - Xl) 2 ) 1/2 

where X is a real solution of the cubic, and Y is the corresponding value from the equa- 
tion for Y given previously. The subroutine finds the real root which gives the smallest 
D. 

For the case of a point solution, the subroutine is specialized with LP = 2, which 
causes the subroutine RESDUE to go to statement 10, and computes 

Y = ( PAZ - AZ(L,K)) 

X = PEL - EL(L,K) 

D = ^X 2 + Y 2 ) 1//2 

where AZ(L,K), EL(L,K) is the data point from site L. 
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Calculation of Root Mean Square of Residuals 

Subroutine SUMRES calls subroutine PAZEL to calculate azimuth and elevation for 
a trial point P, and then calls subroutine RESDUE to calculate the residual from each 
observation site. Then it computes the root-mean-square value of the residuals. 

It usually happens that the cloud viewed at one site does not extend as far as from 
another site. Thus, it is desirable to extrapolate data from such sites, but with reduced 
weighting. Subroutine SUMRES reduces the weighting by the factor 

A 

A + DIST(L) 

where DIST(L) is either zero or equal to DX, the distance between the trial az-el and the 
nearest data point NP(L) on the az-el curve, and A = 0.5/57.3 radians, a somewhat arbi- 
trary fixed angle of 0.5°. Finally, DIST(L) is set equal to zero for the two sites having 
the smallest value of DIST(L). 

Variation of Longitude To Obtain Minimum Residuals 

Subroutine LONMIN(PLAT,PLON,PR,DR,EM) finds the minimum value of the root- 
mean-square residuals EM as the longitude PLON is varied while keeping the geocentric 
latitude PLAT and radius PR constant. 

The procedure is to increment PLON by DA (which is initially equal to DR), chang- 
ing directions when necessary to go through a minimum, calling SUMRES to calculate the 
root- mean- square residual. Consecutive root- mean-square values are labeled (and rela- 
beled as PLON is incremented) El, E2, and E3 so that when E2 is the smallest, the PLON 
which gives minimum can be approximated by using the following analysis. 

Assume that the root- mean- square residual E is a second-order function of the 
longitude DA: 

E = a + bDA + cDA 2 


For minimum E, the longitude is 


DDA = - — 

2c 

From the three values El, E2, and E3, and the corresponding longitudes -DA, 0, 
and DA (using the longitude of the middle point as origin) one can obtain by substitution 
into E 
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El - E3 
b = - ~ 

2DA 

E3 + El - 2E2 
c = 

2DA 2 

from which the longitude (relative to point 2) for minimum root-mean-square residuals is 

(El - E3)DA 

DDA = = — " 

2 (El - 2E2 + E3) 

The increment DA is then decreased by a factor of 10, and the procedure is repeated 
one time. 

Minimization With Latitude 

Subroutine MINISQL(PLAT,PLON,PR,DR,IN,E) starts with the trial point PLAT, 
PLON, PR and while keeping PR fixed, varies PLAT. For each PLAT, LONMIN is 
called to find the minimum root-mean-square residuals with respect to PLON. 

When three consecutive minimum root-mean- square residuals El, E2, and E3 are 
found so that E2 is the least, then an approximation is made (same method as in LONMIN, 
except independent variable is now latitude DR) to determine the PLAT for minimum 
residual. The increment DR is then decreased and the procedure is repeated until IN 
iterations are made (typically three, with consecutive DR values of 0.1°, 0.01°, and 0.001°). 
Upon returning to LARC, the new values of PLAT and PLON give the solution at the par- 
ticular PR. 

Program LARC 

The main program is program LARC, which first reads in the observation site data 
cards in geographic coordinates. For each pass through the program, a triangulation is 
made at one instant of time by using az-el data from simultaneous photographs at each 
site; thus a solution curve of latitude and longitude as functions of altitude is obtained. 

After the time is read in, azimuth and elevation data cards from each site are read 
and stored in AZ(L,N) and EL(L,N); L refers to the site number and N to the data point 
number. Then subroutine SORT(L) sorts the data into either increasing or decreasing 
elevation angles. 
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The BC(I,L,N) coefficients are then calculated for each site L and for every data 
point N by calling the BCOEF subroutine. 

The first estimated solution is chosen near the center of the cloud; since the first 
estimate may be far from the solution, it is better to start in a region where the data have 
the best quality. An integral 100-km value of altitude is used and is read from a data 
card. The first solution point is labeled N = 50 with altitude increments DNA = 100 km 
corresponding to increments of 1 in N, where N now refers to a solution point and not a 
data point. Subsequent trial solutions use the previous solution point. 

One of the main difficulties in the solution occurs near the ends of the cloud. 

Because of differences in exposure, range of cloud, orientation, and visual conditions, the 
cloud visibility may extend farther at one site than at another. Subroutine SUMRES will 
automatically extrapolate curves when necessary, but will give less weight to the extra- 
polated parts. Hence a procedure is needed to stop the calculation when the solution 
curve is going beyond the data from every observation site. This procedure is accom- 
plished with the index IS which is equal to the number of sites using a data point on the 
end of the curve in subroutine RESDUE. When all sites except one are extrapolated 
beyond the end of the az-el curve, then the solution is stopped, and started again at the 
middle (N = 50) with the altitude now incremented downward. 

Solutions are stored in TLAT(N), TLON(N), and TALT(N). For each solution point, 
the point on the az-el curve for each site is stored in NC(L,N) and the root- me an -square 
residuals are stored in E(N) in units of degrees. Finally, the solution is printed out. 

ILLUSTRATIVE CASE 

As an example, actual data from the barium ion release of September 21, 1971, will 
be given. Three observation sites were used: (1) Mt. Hopkins, Arizona; (2) Cerro Morado, 

Chile; and (3) Wallops Island, Virginia. (Coordinates are given at beginning of table H.) 
Table II shows a printout of the az-el input data for the time 3 hrs 18 min 10 sec UT 
(13.307 min after release). There were 36 data points from site 1, 75 from site 2, and 
28 from site 3. The number of points NO used in the least-squares fit was 25, which rep- 
resents about 1.7° of arc when viewed from site 2, and 2.7° when viewed from the other 
two sites. 

Table HI show's the final solution. Total central processor time for the job on the 
Control Data Corporation model 6600 computer was 19.5 seconds. A test was made to 
see how many iterations were made in the solution. For this case, SUMRES was called a 
total of 2922 times, or on the average, 104 times for each solution point at each altitude. 

The root- mean- square residuals were of the order of 0.0045°. 


15 



C ONC HIDING RE MARKS 


A computer program for the triangulation of azimuth- elevation (az-el) data from 
several observation sites has been presented. Because of the relatively short baseline 
used and the high accuracy required, it was necessary to make a special effort to reduce 
the effects of random errors. This reduction was accomplished by using several obser- 
vation sites and many data points from each site. 

An optimal solution was achieved by requiring the minimization of the root-mean- 
square residual of all the sites and by using a least-squares curve fit to the az-el data. 
An illustrative example using three observation sites was given for an actual barium 
cloud. The root -mean -square residuals were of the order of 0.005°. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va. , February 6, 1974. 
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APPENDIX A 


PROGRAM LARC 


Program Listing 

The listing and flow chart for program LARC are presented in this appendix. 


PR OGR AM LAR C l INPUT . OUTPUT I 


* ** 7 HIS PROGRA M SO LVES THE PROBLEM Cf .. L OCATING A CURVED LINE TN SPACE G IVEN 

SIMULTANEOUS AZ-EL DATA FROM SEVERAL OBSERVATION SITES* THE METHOD FINOS 

SJJ CC E S_S 1 V E_ L A T _AN 0 LON POINTS AT SEL ECTED ALTITUDE INCREMENTS. EAC H POI NT 

JS_ VARIED IN BOTH LONG I_T U D E _ AN 0 LATI TUDE IN ORDER TO MINIMIZE TH E RMS 

RESIDUALS FROM ALL THE DOS SITES* _ 

THE RESjnUALS_ ARE COMPUTED US ING A 2ND ORDER L EAST SQUAR ES C URVE FIT TO 

* JHE jiZ-EL DATA IN_THE VICIN ITY O F" EACH TRI AL POINT, 

1 HE PROGRAM El NOS A COMPLETE CURVED LINE SOLUTION FOR SEVERAL EPOCHS 

* OR TJ_w£s _ _ _ __ 

C ”7 ~ _ _J ' ” 

I .< I M E N S |_0 N „ 1'tlAT ( 9 9 ) > T L ON f 9 9 ) , T At T( 99) _ ~ 1“ H'_ 

inMENSJpN , NC I 5 , 99 J _____ 

C UM M 0 N / SI TES/NS , SLAT ( 5J ± SLOM 5 ) ,_SH ( 51 ,UT(5I ,UEG 

_ t GMRf m /l I ME / M ( 5 r 1 9 0 ) , E L ( 5 t I QQ I t BC i 3 , 5 , 1 9 0 ) > NP ( 5 1 » N R ( 5 ) t L P 

L = 6 3 7 i o 2 _ \ P I = 3_o 1^15926535 39S $ RAD=PI/lfl Q t DEG= L/ RAD 

* SET IN ALTITUDE INCREMENT _ _ 

np.iA=ido. 

WT( 1 ) =ViT( 2)=WT(3 I *WT <4 )VWTI 5 >' = !.' _ ' ‘ 

* lead in Dumber cfJh rs sjtes ‘ ' _ _ " " _ ' 

r fid i d,o s 

TERM AT (IS) __ __ 

P Rl NT ~dC)T"" NC "" q'f 6 B s' S I T E S „ " ' “ 1_ 

pri nt I Ql, NS ” _ 

E OR ^ A T <* 1 M 2*“Tf A M ON T R i A N G LJ L A T Y 0 N * / ) 

REA 0 I>r G F 0 G RAP H JC CD" K Q I NA f F S j IN _ DEG REES JU' Q __k vi ) _F n R EACH SI TE AND _ _ 
CON VERT "TO "RADIAN S“~'" ___ 

Dll 11 L = 1_,_NS ' ’ 

L EAD 1 2 f S L A T_j L )_ t_SJLON_( LJ j_S H ( L ) 

Top* a v Xz\- 2o\ i o , V 1 5*2 i 

■RRINJJ 13.2 , SLAT ( L ) f SLCN ( L ) ,SH(L > “ 

F 0 R *"• A T I * SL A T = * F 0* A ^ S LON 0 4 SALT = ,fc FB • 4 ) ’ 

5 L A t TO s SL AT ~(i ) R Ah ~ ~ T S L onTl ) = $ L 0 CM L ) * R A p " 

11 'CPMTINUE_ _ 

******I H lS'j S J3 EG I NNiNG_'OF SOLUTION "FO R EACH E POC H OR T I M E * ~ 

if LEA i) I N T f^f 

1 1 h~ read " 1T1 Yi hF, m i nVsTc 

1 i 1 f Ti P M A T T I I 5 , I 3 1 F 7 • 1 ) 

- j F_ 1 1 f F C ARP l 5 Vl AN K , $T (Vp— T H F” L AS T EPOCH HAS BEEN PROCESSEO*" ' 

iFt I H K o t U * 0 1 S T 'j P _ __ _ 

* ‘PRINT TIME - "" 

PRINT 110 , IHkf _M IN , SF.C 

110 FORMAT I *171 ME* 1 3^ HR^IS* MIN^F^l* SEC^I _ ! ’’ 

* INITIALISE START IMG P n I NT FOR L AT t R S EAR CM DM AZ-EL CURVE 

r p ( 1 1 = NP ( 2 ) = N P ( 3 I =N P ( 4 ) - NP ( 5 ) = 3 

* * * I- FAD I N IJv' P UT M EL DATA F OR EACH SITE 

on 21 


10 

101 

12 

112 
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APPENDIX A - Continued 


* 

r 'f\ D IDENTIFICATION I Nf P R T I ON ON DATA CARO, WHICH HA S COUJM^ 

10 BLANK , 

* 

r N SAME CARD ALSO READ A Z- F L FOR ONE DATA POINT, 



— j 


113 

Ft AD 23fLST, ISl»lS2fTT M EfAZl r Ell 


23 

f OR.* at ( $ x 0 j 46 , A 13 , 22X , 2F 1 1 . 0 > 


* 

IE 1 ST 1 COLUMN 10) IS NOT BLANK, END OF AZ-EL DATA FOR THIS si TF IS' 

±» 

1 NDI CATEP, 



If ( l. S T , N E • j 1 GO TO 2 2 " 



N = M-1 * NblL )=N 


* 

STORE AZ-bL DATA AND PRINT AZ, EL., AMO I DFMT I F 1 CAT I ON INFORMATION 


a Z ( !_ | N i - A Z 1 T ELI L , N J = F L 1 


123 

PR III 123, l, N, AZ( L,NI,FLIL,N1, IS1.IS2, TIME 


L.'.iRVATL* S i A *12* PT * 13 * AZ»#F e .4* EL -* F9.F , A 1 0, A 6 , aToT" 


*: 

Li r\V L RT fcL TO RADIANS ~ 



L L ( L t N ) = 1:. L l L ? N > -PA f) 


jfc 

3'IJT AZ I M PROPER OJA DR AN T AND COVERT TO RADIANS 



IF ( A Z ( L * m J } # 0 T . 1 8 j • ) AZ (L ,N)=AZ(L t N )-36 0* ‘ 



A l (i , r \-,\i ii r rg j - k A •) 



OU T f i U3 ' " 


22 

( ONT I '-I 1 .JP 


*: 

SORT AZ-EL DATA IN p P D t P np pl 


21 

(. AL L. SORT [D 


CUMI^ML ” " 



Ml. U L = 1 , f'iS < NP L = NP I L I 



DO 13 N=1,NBL " 


>\t 

CCMPjrr COFFFICI ENTS FC.R SECOND TRDER LEAST SOUARES'TOCAL Frf 

AT EVERY 

■& 

i‘ATA Pil'iT F'TR EACH OHS SITE. " ' 


1 3 

£ ALL MC^FFILfT ) 


* 

1 L T l l.Ps A IS F ’ R A ! INF SULUMON- IN IS M* OF ITERATIONS IN M INISOL, 

1 Pi. 1 4. T fsi - Ti 


* 1 frAO TIRST f'^Tl '"ATfit PLAT AND' PlrN'T*' ^GGRaPhIc "DEGREES - AT A SELECTED 

* ALII tUOt rt I AND PRINT NOT. ' HI S Hfji.1 L CT ft E~ IN ft 0 HAL MULTIPLES 3 F' A'lT ttude 

* I NI..F £• -1 ENT . — • 

:- LAi) 1 -H » PLA TG » <*L i, n ,HI — 

i:il r ..f "AT (2 Flu. ) ) ’ t 

r p. T ir 11 V, p L ATT, PCT *! , tit 

11a r ''p. v A T ( * L-C ST I V AT F it T»I AL POINT L AT=*F6. 3* loN=*F 7. 3* AO~=*F6 

*<,01 ' . ~ — 

* INITIALIZE '! AND NO Sii rt-AT FIRST S PLOT IP N P 0 \ Mf” I S~ N = _ 5Ci ~A NT) ALT I T UD E HILL 

I OChCASE. 

'■=iv « ND--1 _ “ 

*1. AT I =f’L AT tv* RAD i . . PL ON I = PL ON *~R A [) ~ 

** START IV. PLACE FTP ALTITijDF IMCPFASINC, OR DECREASING _ 

LIS PI AT(>«Pl ATI * °l.r M=PL0NI A N*H I 

*»**,, ^-■.tacT Ing PLACE l-’K EACH SOLUTION POINT 
IF |R=.1«RA0 

* | p £ p rMi:f ; r *■ 

= ^ ♦ gj 

* if f .i v r k t r i r, * n r n t k j c " " 
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APPENDIX A - Continued 


call GCRGCNIPLATG,H,PLAT,PR> _ 

* f I Ml SOL'JTI OM POl NT , I»F. PLAT AND PL UN dHICH MINIMIZES RMS RES IDUAL 

* AT RADIUS PR - 7 ~ 7Z. _ 

.CALL MiNfsoi'i PLAT , PL'i-iN r PR,DR , IM,E< NI_I 

* CP MEBT ' B. AC K Yn GC ngKA PHI C AAIO Of (TRET S AND STIPE ANSWERS . 

.CALL ciCNGGRL PLAT , PR , PL ATG , T ALT I N) » TL AT( M ) =PLAJG*OFG_ _ 

T L " N ( N } - PLC N *0 E G ’ Z"ZZZ~ Z ZZ ... 

« FOR EACH 'SITE STORE OAT A PQ 1 NiT MUM 3 £ P «JH I CH 1$ NEAREST THE SOLUTION 

oc 141 L = 1 »ns " .‘""I...ZZ..Z 

* DETERMINE N'T. t'F SITES WHICH ARE . US l NG._ENO_ .PUL.NI_.DATA .. 

141 •.C(L > M=MPJL ’’ ' _ 

Ts=j ” _ _ 

DO 16 K= 1 1 NS " ’ .. .. 

1 F IN Pit K > . L E 1 1 “ J S •IS ♦ 1 ' " _ _ 

‘ i f i f-i p i k r® g ? » nr i k r> i s = i s + 1 " 

16 CONTINUE ’ 

* LC_ A LL BUT ONE" S I T E 0 S FO E NO - P 0 I NT DATA, JUMP OUT OF LOOP." 

IFUS.GE.NS-1I GO _T0 1 5 _ ... 

* I F PMS RESlOUAL I S T^O~GRTaT JUMP OUT OF LOOP _ 

i fT e i> )“ gt". o7i 3i" go t o ‘ lT" 7_‘ ” __Z Z"_ _ _ _ . 

* IF N IS for LARGE' JEP’p out of loop 

I F ( C , Li). 9 9. CR .au E 0.1 t CO TF iT 7JLZZ1Z 

* 5 RFtiENt" H,'~ GJ TO 14 A : jT ) START VCXT SOLUTI 5N POINT 

f.=H*l)NA*HD " *’ GO TO 14 __ __ 

* IF ALT ITIJOE " I S" HE CREA S j’NG JUMP" ^U T CF "LOOP. ENT IR E LINE $C LUT ION IS 

* NOW FINiSHtO’ " ZZZZ Z Z _ _I 

15 I F I NT) e EOe - 1 > GO T 0 3 31 

* RESET >rTNO_N{r TF START AT CENTER A ND o f c R E A5~F ' a L T I 71 >0 E _ ' 

N2*N "t" ii *51 ~ i N0=- i ’’ ~ _ ' 

^ _ 

331 n i = n 7 7 __ 7 Z ZZZ " ' 

* 'PRINT "time" '"as HEADING FOR COMPLETE SCLUTiON PRINT-OUT 

PRINT 382,1 HRtM IN* SEC" . .... '1. Z1 Z..ZZ _ ... 

382 F0RMATI*1TIMF" «[2^"’hRS " .*12* MIN *#5.1* SEC*) ... 

PRINT 38 . - 

38 f-QRM ATI* LINE SOLUTION*! 

""printTT‘} _ _ „ .... .... 

39 FORMAT (*OAl TIT'i.jr>EZ " LATI T UD E UJNJaJ TIJI) §.. ... RM S RES PTS ON AZ 

3 t-FL CURVE*!' ' '_ _ _ — 

* PRINT OUT TOTAL N O. C-F DA TA PTS, F OR EACH SITE 

PRINT 342 t ( NO ("1)71 =1' . N5 ) ... ... 

3 4 2_ F li RM AT ( 52 X » 5 I 5 I 

V“ ""print OUT SOLUTION POINTS FOP EACH N 

P0_ 34~N = N1 ,N2 ' ^ 

PRINT 36 , T ALT IN) ,~T L AT t N> ,~T LON ( N ) , E ( N > , I NC ( L ,NL , L = 1 ,N S ) 

36 FORMAT I FF » 0 , F 1 2. 3 , F I 3.3, F 1 1 . 4 , 1 1 2 1 4 1 5 I _ 

34 C ONT I NtJE 


GO TO 116 


FMD 
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APPENDIX A - Concluded 
Program LARC Flow Chart 



2 















n : npn,o, 


APPENDIX B 


SUBPROGRAMS 


S U BRUIT I NE. GGRG CN (GLAT.H , CL A T.CR) _ 

*** .SMBPCUTiM. TP CONVERT GEOGRAPH IC LATITUDE GLA T AN D AL TITUDF H 

*** XL'.. ..geoc entric Tat i tuqf c la t. and rau'iuscr 

A =43.7 9 * 1 p 6 3 9 =63 56 . 78 V t E^ B^sTa/ A ~ 

... .. . R= A/S ORT< l. »F* TAN< GLAT 1**21 t_ Z = R* E * T A N ( G CAT ) ^ 

- . *• P*R * H * C [’ S { G L A T I t Zj = l «• HAS I N I GLAT | _ ~~~ 

C R = S Q R T ( S.P.*;R P ♦ Z P *Z P ( __i Cj,AT = AT AN2( zV/RP 1__ ' 

-AETUR^L 1 r. no ~ ~ 


SUbKCUT INE GCNGG R < CL AT, CR tG L AT , H ) 


*** .SUBROUTINE TO CONVERT GEHc'eNT R ic’T AT I t'uDE' CL AT " AND GEOCENTRIC TaTiYuS CR TP 
± G E P GP A P H I C L ATI TUDE G L A T A ND ALTITUDE H. ~ 


. 4 = 63.7 9. 1 06 t F=l,/2«?8.3 ... 

U=CP- Af4*F/2 •*! 1 •—COS C2.*CLAT i + . 5* ( F /ft,- A*F /CR I * IC OS ( A.*CL AT) - l.~) ) 
GLAT = Cl AT *A*F /CR* S I N I 2 *.CL AT i H A*F/CR ) **2*1 1 . -C R/ 4. /A >~*SIN(4.*CLAT ) 

RETURN f. END 


SUB ROUT INC p’azEL { SLAJj SL ON,SH ~. PL AT~. PLON . PR . AZ . EL~I 1.1 ~ 

****GIVJN .AN PBS SITE S ( I N GEOGRAPHIC COORD I . NATFs T AND 

A PL! NI_.P .. < IN GE OCE NTR l.£_ COOR D INATES ) IN SPACE, F I NO THE 

L IN^-OF-S IGHT (AZ,EL>. FROM S TO P. 

A = 63 7 8. 166. . . % .9=6356.784 t.. C=B*B/ A/A " “ 

R»AZSQRf(l.«-E*T AN (SLAT 1**21 i J«R*£*T4.W SLAT L 

* COMPUTE GEOCENTRIC CARTESIAN COMPONENTS. OF. QOS Slif .S 

SR = R*SH*COS_( SLATj. _ % Z S = Z »SU* S I N ( SI. A YY 

XS=SR*COS< SLtiNI * YS = $R*S-I NISLQNI 

* COMPUTE. GEOCENTRIC CARTFSIAN COMPONENTS. L'F POINT P 

XPsPrVc'OS ( PLAT I *COS< PLONI _ _ _ _ 

YP = PR*Ct:s' ( PL AT Y *S IN l PL ON i T ZP = PR*siV(PL AT I 

* COMPOTE COMPONENTS .OF VECTOR. DIRECTI ON FROM SI Y e S TQ~ POINT P 

X A = xjp -X S_ _ _li ~"YA = YP - YS * ZA=ZP-ZS 

* CONVERT TO TOPOCENTR IC CARTESIAN COMPONENTS 

|.X = -X A*SI NISLCN 1 +YA*CPS ISLON ) _ ' 

()Y= - XA* Si M SL 4T ) *CCJS I Si ON) - Y A*S1 N I SLAT ) *S i N ( SLON I F ZA_*COS (SLAT) 
r.Z=*XA*COS(SLATI*CTS (SL.OM)».Y A*Ct>S l SL AT ) *S [Nt SL.CN H ZA *S IN_( St.ATi 

* COMPUTE AZ-EL DIRECTION . J .... 

AZ = A T AN 2 ( l)X j i)Y )_ __ 

Yl=aTan 2 ( Y)'z YsortTox* cx Vdy*dy I ) " _ 

RETURN ' '< END . 1 111 
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APPENDIX B - Continued 


SUFRU'JTINE S'IRT(L) 

C r 1 ■ ■ -1 J 7 L [ N E / A z l 5 , i 9 J ) , £ L ( 5 , 1 90 1 , R C ( 3 ,5 j ISfi A jJjE’C L_ G NB ( 5 I , L P 


I S SUBf OUT INF S.CETS AZ-CL DATA IN EITHER INCREASING OR DECRE ASING VAl UFS 

* if- SI, depending CJM the value of l, so that FIRST value ALWAYS corresponds 

* TO THE HIGH ALTITUDE EMI OF THE CLOUD.. _ 

*** xs j, osT pe modified for a different site configuration hr ''cloud 

* ‘MENTATION. ' ' " “ 


1 

if 


1 


aS*-(-U) 
r r u: i - 1 1 * i 

! F l-L l V A 7 I r r- AM GL 
t XLH4MGF T Mil PME 
J O^CRFASF THt 
s r [J-: j = i » k 

> rtC. ) - N ! i ( L»-l 
T F L - L: L ( L t J ) 

!'L (If U-FL<« ,:j> 

- L ( L j Ml =1 rl 
'•F T s ; l 
■ F L = L L ( L. I J ) 
f. L , J l-fci (L f J + L I 
L ( L * J * 1 ) - T f- 1 
r n.\ T i f « jF 

m-tjfm *, t 


. 1 . M = NB { L J. _ ' 

* K = I " ' 

f r>f n N r point IS EQUAL T Q THAT OF ANOTHER* ’ 
Pni^T .WITH THE VERY LAST pn IMT QN T Hb L I S 
TOTAL NUMBER OF. DATA .P Of MTS hYJo 

* i f a el a r j \ - e tj i.; £ jV i u 2 5 , i cf 

" t " f Ay* ATrrriT 

T AZ(L f J > = AZ<i/N> 

* AZ ( L t *V5 1 = TAZ 


T TAZ-AZ(L ? J) 

= j? z J_Li i±L j 
* AZlL, jn^TAZ* 


*T) 


SUBROUTINE ACOEF(L,N) ... . .. ."’17 

CfJMMON/LINF/ AZ ( 5 , 190 ) , ELI 5 , 190 ) >BC( 3 , 5,1901 , NP (_ 5 > , NBI 5 I ,LP 

l, I ME MS l ON A< 3 , 3 > ,.B <3 »?.C( 50 , 3 J.tipi 3 » 

******* T H I s' SUBMOUTWe "CCMRUTfS T HlTlWrWF FI C 7 E~NTS FOP A SE COND O RD ER 

* LEAST SQUARES CURVE FIT TO THE AZ-EL DATA AT EVERY DATA POINT N. 

* THE C O v PUT AT I CM USES NO DATA POINTS CENTER ED A BOUT N , EXCEPT M THE ENOS 

* V HERE NO END POINTS ARE .USED. THE I NOE PENDENT VA RIAB LE I S ' 

* £ L (l , MV ) -ELI L , N t WHERE fIV IS WITHIN THE S ET ‘is . 

'"'0 = 2 5 “ . . 71 77 

N 0 =N 0»2 ' _ ’ ~ " _ 

to. No=f;o- 2 , .. ".77 

**■ IE THERE ARE FEWER THAN MO DATA POINTS , T HE N DECR EASE N 3 

IF I N B ( L ) • LT « NO ) GO TfVlO " _ __ ‘ ' * 

** ... TH F MIDDLE POINT n~F~ THE SET NO , IS CALLED Nm’~a'N 0 fs' USUALLY' THE SAME AS N» 

«*3 NM=M __ * NC^JNCfU /2 

* AT ENOS OF* CURVE NM I S SPACED NC " PdlNTS" FR 0 ' M ’ THE END. _ ' ’ 

irtH.LT.NC) N**»NC V i Fit N.GT.NB 1 1 ”) Vl-NC ! N M = N B I LI fl-NC 

** THE FOLLOWING COMPUTES THE MATRICES A AND B WHICH ARE USED TO FIND THE 

* LEAST SQUARES FI T . ' " * _ 71 " 

no '20 1 = 1 , NO 

20 ' ; C(T, 1 T=T. _ 

or 5 ) )= 2 ,M i DC SO* 1 = lj NO f> NV=NH -NC»r J ~ 

50 C( I, J l=C( I ,J-II EL( L,NVt-EL(L,N) ) _ _ _ _1 _ 

u n _ f 6 o i = i , m _ * po i oo _ J =;i ,j*_ f_ ~Air t r)=o . I_ ZI 

00 iCC K=i,ND _ _ _ ^ 

L C 0 nj t , J ) = A ( l_i_J UCt K . 1 1 * C < K J | 

up 153 I = 1 ,M ' " % Bi n = 0 . _ 

| ; C 15 0 JS = 1 , \‘f « NV=N M-NC »K 77 Z_ - 

150 MI ) = B('l )tClk,li«( AZ'tL ,N vT- A Z ( L , N II *C 0 S ( ( EM L N ) ♦ EL ( L , N v") j /21 

** filed S'-LVES THF EQUATION AX = B. THE SOLClTiON X ' IS RE TURwifi 'i*N _ fi, WHICH _~‘ 

* CONTAINS THF 3 COEFFICIENTS DESIRED. 

CAL L_s I ±i QI A , M t B ,1 ,D , I P , m, I S ) * ~~ ~ 
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APPENDIX B - Continued 


* 


301 


no 2 ) : i = i , » 

THE CCFFF [C IEKS B AKF STORED IN' BC. 

HC( I * L » N 1= R ( I »” 

! F (M EU.NB(L) r_PRj'NT' 301 * L I'tVijsi 0 

EC R A T T * 0 L = * 12* " ~N=* [ A * ” NQ = * l 3 V 

Ft TURN t gNP " 



' U ;t r. .. i ' f [ - f <E?0*Jr. (!«AZ f PEl.tLi OfXf'J 

;. • 1 V i: j AN El. AND M OF A POINT FROM STATIC 1 N 0, FIND. THE 

*•***■ ; r- r > [Dual ( af distance bfthffu the point and the a;-fl curve) 

i ' N'T./Ll Nu/ AZ ( : ; » i c 0 I i r L t E > 1 *5 0 ) * P» F. | 3 f 5 , IE; I f >'P t 5 .1 , NR (El j£f> . 

C (.•"•o-./SI TLS/^'S, c l AT ( 5 I , SLCNt ; | , SHt 2 I »WT{ 5 | ,DEG 

lF(LP»FJ.3ir,-iTi'ilu 

P '■ IN T f ) HEP U\ CIJKVF CLOSEST TO GIVEN PAZ PEL 

>=; i = ’ :p ( L I T - ! n- 1 

> Si.= ( ( PAZ-AZ ( L f '.‘ I >*C;?S (Df: L ) > ( PEL-ELI L t N> 1**2 .. 

IF |» .GO. MIL I I 00 = -1 

1 .;-Si=nSN . .. 

-m < IFIO.GT.I- rt(L I. OR. I'l. LT.il GO TT. 2 

i SN* 1 I PAZ-AZ U ,M I 1*0 S (Pi III **2* (PEL-EL CLf N) l**i' 

I F ( PS ' L T « PS ” 1 GO To i $ 1FIN.GT.U1- H> GO To 2 

I F ( . L T. Ml GO Z ... 

0 = -l T '>UOO * !>SN = DSP 

T F 1 

2 i -I -Ml * --i P ( L > = M 

****«(-.! VII JTL RESIDUAL FPPN AC COEFFICIENTS CORR E S PGM 01 NG Tfl N 

Al = P£G- r L <L.,NI ...>.. Y 1 = < PAZ-A Z (L ,Nl) l* C iS ( ( PE l+FL I L . N1 I /2 . > 

X 0 = S : J K T ( X 1 n X L * Y 1 s V I J " " ' "" ... 

C ( 1 , L ,N1 T C2 =3C I 2 » L |N I T. C3-B.C1 3,L,N>. 

* i" Zi L G ‘J L A T r THE COEFFICIENTS DF THE CUBIC 

T = 1 . 5*C2 ZC 3 * 0 - < 1. *C2*C2+ ?• *f. 3* ( C 1 — Y 1 ) I Z2.Z C 3/C 3 

• = ( 02 -( Cl-Yl l-Xl I /2./C3/C3 _ 7. .11""..'" 

3. % A = ( 2 . *P*P— F. * Q)* P / 2 7 . FP 

,jC = A<-i'/'2. ♦A*A*A/'2 7." ~ ' 

* IF ill SCR MI»'AFT IS LESS TFIAM ZERO, THERE ARE THREE PEAL SO OTS 

I F t 1C.LE. 0. ) GO TO r ' " 

* •: I SCR (HIM ANT is GREATER THAU ZERO, COMPUTE THE ONF REAL ROOT, 

f DC = Si.it*'T ( i)C1 f CA=-B/2.tRDC * C B = -R_/ 2. -K DC t £=l./3. 

;; =■ S I G r. II A R S I C A 1) * * E , C A I * S |_G Nt ( ABSICF) > **E > CB 1 - PZ3. 

v»ci* CC2*C3 **xT*X i" iis = I Y- Y 1 ) * ( Y-Y 1 I ♦ ( X- X 1 1 * ( X-X11 

•«*SCRUi>S» * PE T'JRN 

f i*rtl=ACOSI-B/’2./S0RT{-4»A*A/27. I I 1 

H S' '= 1 0 0. _ 1 _ ' 

* crp.-PUfE THE THREE REAL ROOTS AND FIUD THE SMALLEST." 

«iC> ! 1 = 1,3 ___ 

/■- 2. s'SORi (-AZ3. i*cV s ['ph!'/ 3'. +f267* U- U /OEGI-PZ 3_. _ _ 

V-C.U (C2*C3*'0*X i PS = ( V- Y 1 V* ( V-Y 1 ) » ( X- X 1 I * 1 X - X 11 

I F (OSO.GT.PS I PSn-DS . _ 

6 UJM T IUUE . - - — 

• l = SOR T (OSC> * RETURN 

****** hip A POINT SJIUMJA LP=2, AMD THE FPLirMNG IS USED 

S = ( PA Z-AZ ( L > K > )*COS( ELI! »K) I ~. . 

A = P L L -E L I t. » k I T D* SORT ( X*X f V*Y I ... 

r F TURN T END . 
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APPENDIX B - Continued 


SUBROUTINE SUHRE S ( PL AT ( PL ON, P R , E ) . 

** JL£P INT » JCAL OJLArrT HV^Rq OT T eAn TsCU AR E OF THE RESIDUALS FROM ~1 

*- "Tt L'' OBSFR V ATl UN SITES. 

*****IE the_point is off the end of t he az-el c urve of a site .then THE nil FACT OR 
'*****''p'^jLL DIMINISH THeVeS [DUAL F ROM THAT SI TE AS DETERMINED BY DISTt THE 

* * * ** A \ G UL A R_0 [ STANCE FROM TflEE NO O F THE C U R V E. HOWEVER, FOR THE TWO SITES 

i****WHlCH. HAVE THE. IE AST ."ui STANCE FROM THE END , make DIS T=Q. SO TH AT F=l, WHIC H 
****wGiyf S FULL WT» _ 

7: ov -iuNTsl tIs /nsTs laTT? > , sloni 5 >» sh Vs i » vit i 5 i t og g ~~ 

COMMON/t INF/ AZI 5,1901 , E L ( 5 > 1 <?0 ) , BC I 3 , 5 , 1 9 0 > ,.NP 15 .) , Nfl I 5 > , L P 

!) INC NS 1 ON 0 1 5 ) > 0 I S t j E . 

a=.f_/oeg r_ - 

nr'i 7=1, ns t jNO-.b. 

* FIND "P AZ'ANI) PE L OF POINT 

• Tali PAZfUSLATI J 1 , SLONI JI ,jtll J 1, PLAT^LON,PR,P_AZ l PEIJ 

* CALCULATE RESIDUAL ' " 

CALL RES but I PA Z, PE L » J , 0 (J V, XD ) .... .... .... 

I F ( UP { J ) . EQ. I .OR » NP ( J ! . EO.NBjJj 1 I " 1 N_0 = 1^_ _ _ _ _ 

L 01 ST I J ) = i Nb«xb 


***** s ; : f<t 01 ST INTO ' 1 NC RE AS I NG 

ORDER 


nft 1Q0 1=1, NS 

i 

K=MS- 1 


l)U I U 0 J = 1 t ^ 


I F t 01 ST ( J ) — 0 1ST ( JUI ) lOOj.100, 1C 

10 

TE«P=biSTC J ) 

* 

DTEM = D< J J 


DISK J > -o I ST C J 1 1 

<; 

DC J»=0( J *1 1 


m ST ( Jfl I =TE”P 

4 

01 Jf 1 ) =DTFM 

100 

CHMT I MiJE 

i 

f = 0 


on 2 L= 1 , NS 

4 

I F ( L * LG* 2 ) ni5T(L)^0* 


F = A /7 AtllTST j' Li r 



<r 

C AL c u late" rm s“r E SlnijAL 



2 

'=E+< OIL) *wf <l ) *F t **2 




E = SDR T’{ E/NS 1 




RETURN t END 




SUoR.VJTlNE Lf OHM PL A tVH6m’,P « , PR 

*.**»**ICK A GIVEN PLAT AND PR, THIS SUBROUJ l.NE F INQS . PLON..WHICH_ GJ.VE.S ,.THE_ 

* "P I Mijm. SMS RE SI iUIAL E v * 

***** -I Nf RF M( NT ol.ON PY i)A LINT I L E2 IS T HE~'l EAST ' OF Th^'t^BE E~ CDMS EC U T I vl 

* H’S RESIDUALS Ei, E2,C3. USING THESE RESIDUALS CO MPUT E A PPROX! MA T E PLPN 

* -Hll.H GIVES THE MI NIMH* (KITH RESPECT TO PLON) RKS RESIDUAL E 1 *, 

* THt;r DECREASE OA And REPEAT THE PROCEDURE, 

!> I*j.l-*1592653589a t 6 E G = lT5 . /7>1 ' ~ ._7"~ 

■ T = t * ft A = D R 

1, : T = E T 1 1 _ 7" I 

■ i - : ' - - 

* FIND K*'S RUSIOU-AL FOP FIRST POIt-'T, 

CALL SiJMocs ( pl at, PL OH, PR, El I J _ 

V I-M1H " _ 

PL. PL'U+UA «. CALL SU^R ES ( PL AT , PLON , PR , E 2 > 

■ I = V I f 1 


1 F U:2 

i * F. 1 F 

(•/: TO 

3 

* 

FS-C.2 £ F 2 = F 1 4- 

F 1-ES % 

FLi M - 

*Li:,- !f ;)A 






PLi L = 

PL Vi + OA 



% 

CALL S U M R E S ( P L AT f PL r 

5 nV PR • F 3 > 

1 ~ ! -i\ 

*1 






IF( S t 

<U.»To2JI 

O r - T-> 

c 




! F ( ■- 3 

CP ) 

G r Tr 

c 

% 

El = : F2 4 F2-E3 « 

GO TO 3 
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APPENDIX B - Concluded 


* FROM T HR ire RESIDUALS FIND Nfc W PL Of* FHR RESIDUAL 

5 . 04s- ( (_ 3 -f 1 1 *04/ t E1-?*F2*-E3 )/?.. * .PL ON* P l ON— 0 A » p D A _ 

lAl'u SUCRES ( PLAT , PLC^f PR> F M >_ _ _ 

IFI'-l .GF.1J)' ' PRINT T»N] ._ 

6 TOP. - - iT(* ITERATIONS IN l DNMI N = * I 4 1 _ 

<,< TP SI fpu 2 ITERATIONS ... . 

]H‘ir,f,p,2i RE torn 

* nt C r L A be $ T F° S ( / [. 

i.-A = D.\/l2. ’ 

i-u r«’ i; 

r f-n 


SUripruTlNF yiNfsrLTPLAT.PLP NtPR. OR. IN, Et . 

***** SUBROUTINE TO PINO A SOLUTION HY .NINiLMI ?ING_ TH.E_pyS_.OF ..THE ,.P_E M mj_AL_$ 

* i'F ALL 'SITE S. W H ILE KEEPING PR FIXED. 

****** 1 k.;c RLMEiJT op AT* BY OP : f f. 4 I 11 NO . 1 ..ON.y IN_ TO F IN D THE MINIMU M R V S OF RE S. 1 DUAL S . 

* (AND THE CORRESPONDING PL 0 N) . _. __ . 

* EC ft" EACH PL AT f UNTIL F 2 IS, TUP | PAST PF THE CONSECU TIVE MINI MUM P^S 

* RF SIGNALS E 1 ,E 2 ,F 3 . ’ " 

* iSir.G THESE RESIDUALS COMPUTE APPROXIMATE PLAT WHICH GIVES THE 


* 

Mi NI 

•MU'I (WITH RESPECT TP DOTH PLAT AND PLONt 

RESIDUAL E. 



* 

t h£T" 

DECREASE. DR AND RC l TFR ATE. 

- - 

- - — - 

•• - 


PI = 3. i 41542653539 3 " E Of 0=130 « /P I 

I T.= •?...... T’_ T OF*. 1. . . _ ... - - 

1 DA = I')I< " i L.I = 0 .... 

* " V INO _F[VST_ P1>N WH ICH MINIMIZES THE RMS RESIDUAL 

call" LCTn^I NC.pl AT , PLONjt PRjOA , £ 1.1... .. _ 

PL 4T= PI AT +0P *. CALL LONMIN! PLAT , PL'TNt PR« D4» E2 I 

IRE2.LT.Fl> GO TO 3 ' * E $ = f 2 F £2*f L. * ..El =3 . S * UR*-”? 

PI. AT = PLAT f DR 

3 f L AT= PL AT ♦OR ’ t CALL L CNM I N ( PL AT , PE ON, P R ,D A , E3 ) 

:. l = L 1_ 

ThTi.ge.25> GO TO 5 ..... . _ 

l FCF3.GT.E2 1 GO TO 5 $ P1=E2 t. L2 = E3..f £G_IQ. 3 

* FRfM_ THREE PESIDUALS FIND NEW PLAT FOR MI NIMUM RMS RESI DUA L. 

5 >DR = - IE 3"E f> * 0 R /' I El - 2 * E 2 ♦ E3 >7f..’ t PL AT = P LAT-D R POOR. 

L A L L L ON N I Nf.PL_AT , P L 6f: * P rT dAjE ) j J I *11*1 * dE«E*D EG ..... 

* TEST F OR NO . CF 1 T ERATI ONS. AND DECREASE STEP. 

i'pljj „EQ, imi Go t q 6 * !y , *OR*OF 1 GO JO _1__ 

6 conTimVie . . - 

>*6e' ’ ... . - ... - 

RETURN « END ... _ 
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APPENDIX C 


COORDINATE SYSTEMS AND CONVERSIONS 

This appendix will describe the coordinate systems used and the subroutines 
involved. 


Earth Model 

For the purpose of triangulating from widely dispersed stations over the Earth's 
surface, the Fischer spheroid was adopted since it is believed to provide the best avail- 
able global fit to the actual geoid. The relevant parameters are 

Equatorial radius; A = 6378.166 km 

Polar radius: B = 6356.784 km 

A - B j 

Flattening factor: F = 

~ A 298.3 


For this model, the deflection of the vertical, that is, the angle between the normal 
to the geoid and the normal to the Fischer spheroid nowhere exceeds 30 arc seconds which 
is sufficiently accurate for the present purposes, since pointing directions used in trian- 
gulation are referenced to the stars rather than to a local horizon. 

Geographical Coordinate System 

The geographical coordinate system is the conventional system of latitude, longitude, 
and altitude. Figure 6 shows an exaggerated spheroidal surface corresponding to the 
Earth’s sea-level surface. The geographic latitude GLAT of a point P is the angle between 
the equatorial plane and a line drawn from P perpendicular to the spheroidal surface. The 
altitude H is measured from the surface at point G to the point P. The longitude GLON is 
measured eastward from Greenwich. 

Geocentric Coordinate System 

In this system, O is the Earth's center in figure 6. The X-axis is directed toward 
the intersection of Greenwich meridian with the equator. The Y-axis is directed toward 
90° east longitude in the equatorial plane, and Z is directed toward the north geographic 
pole. The point P is also located by the geocentric longitude GLON, geocentric latitude 
CLAT, and radius CR from the Earth’s center. 
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APPENDIX C - Continued 


Topocentric Coordinates 

This is a local system with center at some observation site S (see fig. 4) with the 
DX,DY plane coincident with the horizontal plane, with DX directed toward east, DY 
directed toward north, and DZ directed vertically (that is, perpendicular to the surface of 
the spheroid). In the polar version, a point P is located by azimuth angle AZ measured 
clockwise from DY (north), and elevation angle EL measured up from the horizontal plane, 
and range RA measured from S to P. 

Conversion From Geographic to Geocentric 

This conversion is accomplished by using subroutine GGRGCN. Reference to 
figure 4 shows that the point G on the spheroid follows the equation 

R 2 Z 2 



where 

R 2 = X 2 + Y 2 


The slope on the ellipse is 


dZ 

dR 


B 2 R 


A 2 Z 


hence 


tan GLAT = 


B 2 R 


or 


Z = R E tan GLAT 


where 
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APPENDIX C - Concluded 
Substitution in the original equation gives 

A 

5 = 

|l + E tan 2 GLAT 

The geocentric coordinates of point P are then 
RP = R + H cos GLAT 
ZP = Z + H sin GLAT 


CR = 

CLAT = tan-! = 

CR 

Geocentric to Geographic Conversion 

This inverse conversion cannot be obtained explicitly. Subroutine GCNGGR uses 
the following approximate formulas derived in reference 12. 

The altitude is given by 

J ,/f a F\ / \ 

H = CR -A + -AF1- cos (2CLAT) + cos (4 CLAT) - 1 

- — - 2 2\4 CR J\ 1 

. 

where CLAT is the geocentric latitude and CR is the geocentric radius. 

The geographic latitude is given by 



A F 

GLAT = CLAT + — - sin (2CLAT) + 
CR 



cr\ 

~|sin (4CLAT) 
4A 


It should be noted that all angles must be expressed in radians. 
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TABLE L- COMPARISON OF VARIOUS TRIANGULATION METHODS 


Reference 

Author’s 

name 

Type of 
object 

Number of 
observation sites 

Method of solution 

2 

Brown 

Point 

Many 

Least- squares method 

3 

Hogge 

Point 

Straight line 
Curved line 

Several 

2 

2 

Least- squares method 
Intersection of ray from one site with 
plane from other site 
Intersection of ray from one site with 
surface defined by 3d order least- 
squares fit to data from other site 

4 

i 

i 

Lloyd 

Points 

Straight line 
Curved line 

2 1 

2 

2 1 

! 

Midpoint of minimum skew distance 
between rays from each site 
Intersection of two planes 
Intersection of ray from one site with 
surface from other site 

5 

Whipple 

Straight line 

2 

Intersection of two planes 

6 

Justus 

Point 

Curved line 

2 

2 

Equal residuals 

Intersection of ray from one site with 
surface from other site 
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TABLE II. - INPUT DATA FOR ILLUSTRATIVE EXAMPLE 


3 STATION TR I AI GUL AT ION 


SLAT a 

■31. 

68 5 3 

SLON=- 11 0.9774 

SALT = 2.3640 



SL AT = 

-3G. 

1666 

SLON= -70.7673 

SALT= 2.1346 

• 

SLAT- 

37c 

9324 

SlGN = -75.4717 

SA 

LT= ,0106 

' 

time 


3 HR 18 

MIN 10.0 SEC 





ST A 

1 

PT 

1 

A 2= 119.0452 

E L- 

43.2735 

MT HCPKINS CA 1 

0-2 3,13,1 

STA 

1 

PT 

2 

A Z = 119*1804 

EL - 

43.2085 

MT HOPKINS CAM 

C-2 3*18,1 

STA 

l 

PT 

3 

AZ= 119.2913 

EL = 

43.1431 

MT HOPKINS CAM 

C-2 3,18,1 

STA 

1 

PT 

4 

AZ = 116.4330 

EL - 

43.0773 

MT HOPKINS CAM 

C-2 3.18.1 

STA 

_x_ 

PT 

5 

AZ = 119.5394 

EL- 

43.0121 

MT HOPKINS CAM 

C-2 3,18,1 

_ STA _ 

l 

? T 

6 

AZ = 119.6643 

EL* 

_.4.2 .9467 

MT HOPKINS CAM 

c_t£ 3_,_ie,_i. 

STA 

1 

PT 

7 

AZ = 119.8046 

EL - 

42.3324 

MT HOPKINS CA M 

C-2 3 jl 1 8_t I 

STA 

1 

PT 

8 

AZ = 119.9364 

EL- 

42.3175 

MT HOPKINS CAM 

C-2 3^13,1 

STA 

i 

PT 

r 

AZ = 120.070? 

FL = 

42.7521 

MT HOPKINS CAM 

C-2 3,18,1 

ST A 

l 

PT 

id 

A Z = 120.1341 

EL- 

42.6963 

M T HOPKINS C A? 1 "- 

C-2 3 ♦ 18*1 

STA 

i 

PT 

11 

A Z • 1 2 C * 3 1 ID 

EL- 

42.6210 

MT HOPKINS CAM 

C-2 3,18,1 

STA 

i 

PT 

12 

AZ* 120.4449 

EL- 

42.5552 

MT HOPKINS CAM 

C-2 3,18,1 

STA 

i 

PT 

13 

A7 = 120.5667 

EL - 

42.4993 

MT HOPKINS CAM 

C-2 3,18,1 

STA 

x 

PT 

14 

AZ= 120.678o 

EL- 

42.4229 

MT HOPKINS CAM 

C-2 3,18,1 

STA 

l 

PT 

IV 

AZ - 120. 7976 

FL = 

42.3563 

MT HOPKINS CA* 

C-2 3,19,1 

sta 

l 

PT 

It 

AZ = 1 20. c 22 3 

EL* 

42.2399 

MT HOPKINS CAM 

C-2 3, 19,1 

STA 

l 

PT 

17 

AZ = 121.0436 

FL- 

42.2229 

MT HCPKINS CAM 

C-2 3,18,1 

sta 

l 

PT 

13 

AZ = 121. 1 7 1 4 

EL- 

42.1565 

MT HOPKINS CAM 

C-2 3,13,1 

STA 

l 

PT 

19 

42= 121.2312 

EL- 

42.0399 

MT HOPKINS CAM 

C-2 3,13,1 

STA 

i 

PT 

20 

AZ - 121.4022 

FL- 

42.0226 

MT HOPKINS CAM 

C-2 3,18,1 

STA 

l 

pr 

21 

AZ' 121.5237 

EL- 

41.9556 

MT HOPKINS CAM 

C-2 3,13,1 

STA 

i 

PT 

22 

A Z ' 121.6460 

FL- 

41.8837 

MT HOPKINS CAM 

C-2 3,19,1 

STA 

i 

PT 

23 

AZ = 121.7643 

EL- 

41.3212 

MT HOPKINS CAM 

C-2 3,18,1 

sta 

i 

PT 

24 

AZ- 121.8690 

EL- 

41.7536 

MT HOPKINS CAM 

C-2 3,18,1 

STA 

i 

PT 

25 

AZ- 121.0944 

EL = 

41.6063 

MT HOPKINS CAM 

C-2 3, 18,1 

STA 

i 

PT 

26 

AZ - 122.1165 

Ft- 

41.6182 

MT HOPKINS CAM 

C-2 3, 18,1 

STA 

i 

PT 

27 

A Z - " I 2 2 m 2 2 7 1 

EL- 

41.5496 

«T HOPKINS CAM 

C-2 3,18,1 

STA 

i 

PT 

28 

AZ- 1 22 o 3316 

FL- 

41.4809 

MT HOPKINS CAM 

'C-2 3,18,1 

STA 

i 

PT 

29 

AZ = 122.462? 

EL- 

41.4121 

MT HOPKINS CAM 

C-2 3,18,1 

STA 

i 

PT 

30 

. - - 

AZ = 122.5736 

EL- 

41.3435 

MT HOPKINS CAM 

C-2 3,18,1 

STA 

i 

PT 


AZ = 122.6364 

FL- ' 

“41.2 748 

MT HOPKINS CAM 

C-2 3,18,1 

STA 

i 

PT 

"if 

' AZ = T 22 » fTo 7 4 

EL- 

4 1. 2 057 

MT HOPKINS CAM 

c-2 3, i3, r 

STA 

l 

PT 

..... 

AZ - 122.9U5 

FL- " 

41". 13 68 

M T HOPKINS CAM 

C-2 3,19,1 

STA 

i 

PT 

34 

A 1 = 123. 0 3 7 9 

EL - " 

4 1.0 674 

MT HOPKINS CAM 

C-2 3,13,1 

STA 

i 

pf 

35" 

■ A Z = 123. 1260 

FL- 

40.9998 

MT HOPKINS CAM 

C-2 3, 18,1 

STA 

i 

PT 

3 6 

A Z = 123.2700 

FL- 

40.9294 

Mf HOPKINS CAM 

C-2 3, 18, 1 

STA 

2 

PT 

"’1X7 

AZ - 3 49.5 33 5 

FL- 

'“49. 12 04 

CHILE CAM ~D-3 

3,18,10 

STA 

2 

' Pf 

2 

A Z = * 3 49 « 5 51 9 

FL- 

4 9 0 C 5 2 4 

CHILE' CA« D-3 

3,18,10 

S TA " 

2 

PT 

3 

AZ = 3 49«5:>2 5 

F L - ~ 

48.9852 

CHILE' CAM 0-3 “ 

3,18,10 

STA 

2 

PT 

4 

A 7= 349.5 613 

EL - " 

48.9175 

CHILE CAM' D— 3 

3 ? 1 8,10 

STA 

2 

Pf‘ 

c 

AZ = 349.5991 

EL- 

49.9493 

CHILE CAM 0-3 

3,19,10 

STA 

2 

Pt 

t " 

A Z =' 34 9. 535 5 ‘ ' 

FL- 

4F. 7822 

CHILE CAM 0-3 

3 , 19 , 1 0 " 

STA 

2 

PT 

7 

A Z = 34 9.616 0 

PL.- 

4 9.7143 

CHILE CAM 0-3 

3,18,10 

STA 

2 

PT 

3" 

AZ- 349.6280 

EL- 

‘49.6468“ 

CHILE CAM 0-3 

3,18,10 

STA 

2 

PT 

9 

AZ- 3 49.5405'“ 

E L - 

'4 9. '5 7 39 

CHILE CAM D- 3 

3,15,10 

STA 

2 

PT 

10 

AZ- 349.635 0 

FL - ' 

49.5114 

CM I L ECfAM 0-3 

3,19,10 

STA 

2 

PT 

"if 

AZ- 34 9.6699 

EL - " 

48.4436 

CHILE CAM 0-3 

3,18,10 

STA 

2 

PT 

12 

A 2-* 34 9. 685 9 

FLil 

40.3762 

CHILE CAM 0-3 

3,18,10 

sta 

2 

PT 

13 

A Z -"34 9.6937 

EL- ’ 

48.3075 

CHILE CAM 0-3 

3,19,10 

STA 

2 

PT 

14 

AZ = 34^. 7015 

EL - 

48.2403 

CHILE CAM D— 3 

3,19,10 

STA 

2 

PT 

15 

AZ- 349.7223 

FL- 

48.1732 

CHILE CAM 0-3 

3, .18, 10 

STA 

2 

PT 

16 

AZ= 349.7396 

EL- 

48.1047 

CHILE CAM 0-3” 

3,18,10 

STA 

2 

PT 

17 

A Z “ 3 +9.7 54- 

CL- 

48.0372 

chile cam 0-3 

3,18,10 . 

STA 

2 

p T 

18 

AZ = 349.7631 

FL- 

47.9686 

CHILE CAM 0-3 

3,18,10 
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TABLE II.- INPUT DATA FOR ILLUSTRATIVE EXAMPLE - Continued 


STA 

2 

PT 

19 

AZ = 3 4'i„7 79? 

Ft * 

47.9017 

CHI LF 

CAM 

0-3 

3, 18,10 

jTA 

2 

PT 

2 

AZ* 349.791ft 

EL* 

47.8341 

CHILE 

CAM 

D— 3 

8,18,10 

STA 

2 

PT 

21 

AZ = 349.3059 

Ft * 

47.7669 

CHI LE 

CAM 


3,13,10 

STA 

2 

PT 

? 2 

A 7 = 360,3234 

EL * 

47.6998 

CHILE 

CAM 

0-3 

3, 18 ,10 

STA 

2 

PI 

23 

A Z = 349.3379 

Ft = 

47.6326 

CHILF 

CAM 

0-3 

3,18,10 

STA 

2 

PT 

2 A 

AZ* 349.8514 

Fl* 

47.5647 

CHILE 

CAM 

D— 3 

3,18,10 

STA 

2 

PT 

25 

A Z - 34°. 962ft 

EL* 

47.4974 

CHILF 

CAM 

0-3 

3,18,10 

STA 

2 

PT 

2 6 

A Z = 349**782 

H L = 

47.4299 

CHILE 

CAM 

D— 3 

3, Ifl ,10 

STA 

2 

PT 

27 

A Z = 349.5911 

EL* 

47.3620 

CHILE 

CAM 

0-3 

3,18,10 

STA 

2 

PT 

? p 

A Z “ 349.9025 

Ft - 

4 7.2 944 

CHILF 

CAM 

n-3 

3,13,10 

STA 

2 

PT 

2 r 

A L- 349.9153 

FL* 

47.2261 

CHILE 

CAM 

!>— 3 

3, 19 ,10 

STA 

2 

PT 

3 5 

A Z = 3 49.92 86 

Ft - 

4 7 0 1 5 09 

CHILF 

CAM 

0-3 

3, 18,10 

STA 

2 

PT 

31 

A Z = 349.9413 

EL' 

47*0913 

CHILE 

CAM 

0-3 

3, 19,10 

STA 

2 

^ T 

32 

M- 340.9 535 

EL = 

47.0237 

CHILF 

CAM 

D-3 

3,18,10 

STA 

2 

PT 

33 

A Z = 349.9677 

EL* 

46.9561 

CHILE 

CAM 

0-3 

3, 18, ID 

STA 

2 

PT 

34 

A Z = 34 9. ''9 0ft 

EL * 

46.9992 

CHILE 

CAM 

n-3 

3,18,10 

STA 

2 

PT 

35 

A Z - 345.5954 

EL = 

46.8210 

CHI LE 

CAM 

0-3 

3,18,10 

:> TA 

2 

PT 

2 0 

AZ = 350.013° 

EL' 

46.7535 

CHILE 

CAM 

0-3 

3,18,10 

STA 

2 

PT 


A Z = 3' Jo j 245 

FL = 

46.6355 

CHILE 

CAM 

D-3 

3,19,10 

STA 

2 

PT 

33 

A/ = 3 5 0. )3 30‘ 

Fl* 

46.6i82 

CHILE 

CAM 

0-3 

3, 18, 10 

STA 

2 

PT 

3 9 

A Z = 

EL = 

4 6 0 5 5 06 

CHILE 

CAM 

0-3 

3,19,10 

STA 

c 

PT 

4J 

A Z = 35D,D52c 

EL - 

46.4829 

CHILE 

CAM 

D-? 

3,19,10 

STA 

2 

PT 

41 

A Z = 350. Jo 5 5 

C L = 

4 t . 4 1 46 

CHILF 

CAM 

n-3 

3,18,10 

STA 

2 

PT 

42 

AZ = 33 0. 3785 

EL* 

46.9474 

CHILE 

CAM 

D-3 

3,19,10 

STA 

2 

PT 

43 

4 7= 353. }5C3 

FL = 

46.2798 

C H I L F 

CAM 

0-3 

3,18,10 

STA 

2 

PT 

*T 4 

A Z - 33j.1o99 

EL = 

4o.2123 

CHILE 

CAM 

D-3 

3, 18, 10 

STA 

2 

PT 

45 

A Z = 350.1133 

EL* 

46.1446 

CHILE 

CAM 

n-3 

3,19,10 

STA 

2 

PT 

* f 

AZ= 350.1340 

EL = 

46.0768 

CHI LF 

CAM 

D-3 

3,18,10 

STA 

2 

PT 

47 

A Z “ 350.141ft 

FL * 

46.0091 

CHILE 

CAM 

D-3 

3, 18,10 

STA 

2 

PT 

4 4 

A7* 3fO.l506 

EL' 

^5.9418 

CHILE 

CAM 

0-3 

3,18,10 

STA 

2 

PT 


A Z r 350.1643 

FL* 

45.8742 

CHILE 

CAM 

0-3 

3,19,13 

STA 

2 

PT 

5 u 

£7 = 35 3.1784 

EL* 

4 5 « 8 £7 3 

CHI LE 

CAM 

0-3 

3,18,10 

STA 

2 

PT 

51 

aZ= 353.1QG5 

EL - 

4 5. 7 3 84 

CHILE 

CAM 

0-3 

"3,18,10 

STA 

2 

PT 

52 

AZ = 3 5 Jo 2 « 3b 

f t = 

45.6705 

chile’ 

C AM 

0-3™ 

~ 3,18,10 " 

STA 

2 

PT 

C. T 

A Z ~ 3 r 0 □ 2 1 5 6 

r L * 

4 5.6 02 5 

CHILE 

CAM 

o-T 

”3,18,10 

STA 

2 

PT 

54 

AZ= 3 5 C d 2 2 7 2 

F l - 

4 5 . 5343 

CHILE 

CAM 

“0-3 

3,18,10 

STA 

-> 

L. 

PT 

5 r 

AZ= 3<0.2 333' 

rc = 

45.4663 

CHILE 

CAM 

n-I' 

3 , IB , i 0 ' 

STA 

2 

PT 

’ 56 ‘ 

’ 47= T30;TTj2 ’ 

rev' 

45.2994 

CH I LF 

CAM 

D-3 

3, 19,10 

STA 

2 

PT 

57 

AZ = 3 rj. 27-37 

EL = 

4 5'. 3 30T 

’ CHILE' 

Tam ’ 

D-3 

3, 18,10 

SI A 

2 

PT 

5 ft 

AZ = 3 50.2 761 

cL = 

45.2620 

CHILF 

CAM 

n-3 

3,19,10 

STA 

2 

PT 

5; <0 

A 7. - 3' 0.2 370 

EL*' 

' 4 5 , 1924 

CHILE 

CAM 

D-3 

3,13,10 

sta 

2 

PT 

tv 

A Z= 3 5 ^ .2° 93' 

FL - 

4 5. 1255 

CHILE" 

“cam" 

0-3 

3,18,10 

STA 

2 

PT 

61 

\2 r 35 0;304T 

F L = 

"■ 4T.05T2 

“CHILE 

CAM” 

n-3 

3, 18 , 10 

STA 

2 

PT 

O c 

AZ= ’35 0. 3 196 

FL* 

4^.0396 

CHILL 

CAM 

0-3 

3, IS, 10 

STA 

2 

PT 

t-3 

A Z *" 3 5 0.330° 

FL*” 

44.9218 

CHILE" 

CAM 

0-3 

3,18,10 

STA 

2 

PT 

6 4 

A l - 3 5 J . 3 4 b 5 

r. L - 

44.9539 

'CHILE 

"C AM 

0-3 

3,19,10 

STA 

2 

PT 

65 

AZ = 35 0.3578 

FL - 

44.7854 

chile 

CAM 

D-3 

3,18,10 

STA 

_2„ 

PT 

66 

A Z =__3_5. 0 . 3 6 99 

FL_r . 

44.7175 

CHILF 

CAM 

_DrJ. 

3,18.10 

sta" 

2 

PT 

.67 

A Z = 350.3859 

FL=L 

44.6494 

CHILE 

CAM_ 

0-3 

3,18.10 

STA 

2 

PT 

6 r> 

AZ= 35O.395I 

FL- 

. 44.8911 

CHILE 

C A M 

D-3 

3,19,10 

’ STA 

2 

pT 

”69 

A Z = 353.4368 

FL* 

44,5129_ 

.CHIT E 

' cam" 

D-3 

3,18,10 

STA 


" pt' 

1.0... 

. AZ= 35 0.42.07 

EL* 

44.4445 

CHILE 

CAM 

0-3 

3,13.10 

STA 

2 

PT 

71 

AZ =_ 3 5 J. 43 13 

FL- 

44.3759 

chile 

CAM 

0-3 

3,19,10 

STA 

2 

PT 

72 

AZ =■ 350.4443 

EL - 

'*4.3371 

CHILE 

CAM 

D-3 

3 , 1 9 ,_1 0 _ 

ST A ' 

2 

PT 

73 

AZ = 350.4586 

FL- 

44.2390 

CHILE 

CAM 

D-3 

3,13,10 

STA 

2 

PT’ 

74 

A Z = 3 50. 4677 

FL - 

. ..44.17 04 

CHILE 

j5AM_ 

0-3 

3, 19,10 

“sta 

2 

PT 

75 

A Z - 350.4737 

FU- 

44.1016 

chile” 

CAM 

0-3 

3,13,10 

sta 

3 

PT 

1 

A Z = 1B2.3343 

EL* 

50.5829 

AC WALLOPS 

i AC- 

-3 3,18.10 

STA 

3 

PT™ 

2 

AZ = 13 2.3085 

EL - 

50.6999 

AC WALLOPS 

, AC- 

-3 3,1S, 10 

STA 

3 

PT 

3 

AZ= 182.3122 

"fi- 

50.7940 

AC WALLOPS 

i - AC- 

-3 3,18,10 

STA 

3 

PT' 

4 

A Z = 182.3192 

ll- 

50.0019 

AC WALLOPS 

r'AC" 

*3 3,19,10 



TABLE n.- INPUT FOR ILLUSTRATIVE EXAMPLE - Concluded 


STA 

3 

PT 

' R 

AZ = 18 2. 3 341 

el- 

51.0146 

AC 

WALLOPS 

AC -3 

3,18,10 

ST A 

3 

PT 

6 ' 

AZ = 182.3490 

EL - 

51.1274 

AC 

WALLOPS 

AC-3 

3,18,10 

sta 

3 

PT 

7~ 

AZ = 182.3511 

FL = 

51.2313 

AC 

WALLOPS 

AC- 3 

3,19,10 

STA 

3 

PT 

A 

AZ - 182.3554 

EL = 

51 .3366 

AC 

WALLOPS 

AC-3 

3,18, 10 

STA 

3 

PT 

a 

A Z = 132.3601 

‘ El = 

51.4417 

"AC 

WALLOPS 

AC -3 

3,18, 10 

STA 

3 

PT ' 

10 

' AZ ="192,3723 

EL ” 

51.5526 

AC 

WALLOPS 

AC-3 

3 , 1 fl , 10 

STA 

3 

PT 

11 

A Z = 132.3321 

FI = 

51.6622 

AC' 

WALLOPS 

AC-3 

3,19, 10 

sta 

3 

PT ' 

12 

AZ " 132. 3391 

EL = 

' 51.7683 

AC 

WALLOPS 

AC-3 

3,18, 10 

STA 

3 

PT 

13 

" AZ = 132,3 95 3 

‘fl = 

5 1 o 3 749 

' AC 

WALLOPS 

AC-3 

3,18, 10 

STA 

3 

PT 

14 

A Z = 182.4013 

F L = 

51. 9808 

AC 

WALLOPS 

AC-3 

3,13, 10 

STA 

3 

PT 

15 

AZ^ 132.4039 

EL = 

52.0844 

AC 

wallops” 

AC-3 

3,18, 10 . 

STA 

3 

PT 

16 

AZ = 132.4213 

EL = 

52. 1985 

AC 

WALLOPS 

AC-3 

3,18, 10 

sta 

3 

PT 

17 

A Z= 192.4277 

EL = 

52.3046 

AC 

WALLOPS 

AC-3 

3,18,10 

STA 

3 

PT 

14 

AZ = 182,4363 

EL = 

52.4123 

AC 

WALLOPS 

AC-3 

3,13,10 

STA 

3 

PT 

19 

AZ = 152.4439 

EL = 

52.5196 

AC 

WALLOPS 

AC-3 

3,13,10 

sta 

3 

PT 

20 " 

4Z= 182*4526 

EL- 

52.6272 

AC 

WALLOPS 

AC-3 

3,18,10 

sta 

3 

PT 

21 “ 

AZ= 132,4655 

EL - 

52.7376 

AC 

WALLOPS 

AC-3 

3,18, 10 

STA 

3 

' PT 

2 2 

AZ = “1 R2.‘4‘7 09 ' 

' EL - 

"52.8434 

"AC 

WALLOPS 

AC -3 

3,18,10 

STA 

3 

PT " 

2 3 

AZ = 132.4764 

'"' EL* 

52.9492 

""AC 

'WALLOPS' 

AC-3 ' 

3,18,10 

STA 

3 

Pt 

24 

AZ= 132.43 72" 

EL - 

" 53.0 530 

"AC 

"wallops 

AC-3 

3 , i'8,10 

STA 

3 

PT 

? 3 J 

' AZ = 1 82.4923" 

L L = 

53* 1 636 

AC 

WALLOPS 

' AC -3' 

"3,18,10 

STA 

3 

PT 

“It" 

AZ = 1*2.5001 

EL- 

53.2 70S 

AC 

WALLOPS 

AC-3 

3 , 1 8 i 1 0 

STA 

3 

PT 

2 7 

A Z = IH2 .5 079“ 

EL = 

53.3 786 

" AC 

WALLOPS 

AC -3 

3,1 8,1 0 

STA 

3 

PT 

28 

A Z = 132.5226 

" ~CL * 

53.4909 

AC_ 

'wallops 

AC-3 

3,18,10 













L= 1 


N = 

'36 ' 

N0= 25 


















L =s 2 


;>! ~ 

7 5 

_M 0 *_ J.5 


















L = 3 



26 

"NO = ^5 " 


















estimated trial 

POINT L AT - 6 . 

15C_” 

rON^-76". 

770' 

;"altt_' 

31500 - 
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TABLE III.- SOLUTION FOR ILLUSTRATIVE EXAMPLE 


t I h' t; * HRS 

13 V-I\i 

1 v • 0 5 5 C 





L INI SOLUTION 



■ • 

- 


- - 

ALTITUDE l 

-m tu ’)*: 

lONCA TiJDc 

RMS PES 

PTS ON" AZ- 

E L 

curve” 





36 

75 

28 


. l ? i 

-76* 7 C 3 

. 0024 

35 

75 

23 

3030 0 

7 » 9 7 ^ 

- 76*763 

.0024 

34 

72 

23 

30 4L ^ 

7 m 83 6 

-7*. 763 

.002 5 

32 

69 

28 

30 3uo 

7 o t f ■ o 

- 76* 762 

.0027 

31 

£7 

28 

306u - 

* • 55 6 

- 7 6. 76 7 

.0027 

3C 

64 

28 

3070 0 

7 .416 

-76. 762 

.0033 

29 

61 

2 8 

3 j 8 0 j 

/ A <L. { 9 

- 76 # 7 6 7 

. 0039 

27 

5 Q 

23 

3 09 J 0 

7*137 

-76. 76 1 

.0040 

26 

56 

2 7 


7. jj3 

-76. 0 

. 0044 

25 

r >4 

25 

3UU0 

7.369 

-76. 758 

.0046 

24 

51 

2 4 

312 0 0 

i o 7 2 8 

-76. "55 

. JC44 

22 

43 

23 

3 1 30 

< • 6 v L 

-76.75 6 

.0 04 7 

21 

46 

22 

3u: ; 


- 77; 75 5* 

.0051 

20 

44 

20 

315)0 

D J C 0 

-76.754 

. 00 54 

" 1 9 " 

41 

19 

3160 3 

• 13 2 

-76.75 2' ' 

.0057 

17 

"3R 

' """18 

J 1 70 j 

C • v* J 9 

-76. 751 ' 

• 00 5 B 

it 

35 

16 

3130 0 

■so ! s 6 O 

/ C o i z> 0 

. l C 6 1 

“15 

3 3 

15 

5190 j 

Vo 751 

-76.749 

. 006 2"' 

13 

30 

13 

32 J J J 

:■ . 0 3 

""- tT. 7 4 5 " 

o 00 5 8 

12 

27 

12 

J.21U0 

5.44* 

-76.747 

• 0055 

1“ 

24 

!! 

3220 0 

: * 310 

-7t.746 

o 00 5 3 

9 " 

' 22 

9 

32300 

• . 1 1 1 

-7 6.74 5 

"V0D4 8 

8 

19 

3 

32 4u J 

' .011 

-7 6. 74 6 

"*^042" 

”"T" " 

'16 

6 

32 500 

- o u t 7 

-76 . 744 

• 0042 

K 

' i4 

5 

3260 -j 

~o 775 

-76 7f?2 

. 0 U 4 h 

4 

11 

4 

327C J 

' m 034 

-76.741 

. 004 8 

‘ 2 

8 

2 

J>2 3„ j 

4.4^7 

- 7 6* 7 3 9 * 

© 00 52 

T 

... & 

1 
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PLAT 

C 
B 
A 

1 PLON 

Figure 3.- Illustration of minimization of root- mean-square residuals with 
varying PLON and PLAT. EM is the minimum root- mean- square 
residual for variable PLON at a fixed PLAT. 
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(NORTH POLE) 
Z 



(GREENWICH) 

Figure 6.- Illustration of the Earth spheroid and the relations between 
geographic and geocentric coordinates. 


40 


NASA-Langley, 1974 


L-9268 


